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BACKGROUND  OF  THE  JCF  RULE 


Bak,  Tang,  and  Weisenfeld  (BTW)  (ref  1)  introduced  the  concept  of  self-organized 
critically  (SOC)  to  provide  a  consistent  explanation  for  the  fractal  spatial  structures,  power-law 
distributions,  and  flicker  noise  commonly  observed  in  spatially  extended,  dissipative,  dynamical 
systems.  Jensen,  Christensen,  and  Fogedby  (JCF)  (ref  2)  clarified  the  ideas  in  BTW  and 
established  the  connection  between  the  power-law  dependencies  of  the  distribution  of  lifetimes 
and  the  power  spectral  densities  (PSD). 

Denoting  the  lifetime  of  an  event  by  T,  the  JCF  weighted  distribution  of  lifetimes  G(T)  is 
defined  as 


GO)  =  I  dSP{S. 

0 

where  P(S,  T)  is  the  joint  probability  for  total  time  integrated  "sliding"  S  and  lifetime  T  of  an 
event  (e.g..  an  avalanche).  JCF  demonstrate  that  the  PSD  corresponding  to  such  a  weighted 
distribution  of  lifetimes  is  given  by 


Sif) 


— ^  f  drG(r)smHnJr) 


(2) 


where  v  is  the  pulse  repetition  rate  and  /  is  the  frequency. 


JCF  assumed  that  the  weighted  distribution  of  lifetimes  C(T)  varies  (approximately)  as 


G(7) 


fr-exp(-77ro). 


when  T  <  Iq 
whenT  2. 


(3) 


We  refer  to  the  form  of  Eq.  (3)  as  an  "exponentially  cutoff^  power-law  distribution  of  lifetimes." 
The  parameters  in  the  distribution  are  referred  to  as:  (1)  the  lifetime  distribution  exponent:  a; 
(2)  the  exponential  cutoff  parameter:  T^;  and  (3)  the  lower  lifetime  cutoff:  V 

JCF  established  that  an  exponentially  cutoff  power-law  distribution  of  lifetiMcs  gives  rise 
to  PSD  S(/)  of  the  form. 


S(/)  « 


const. 


when  f  <  l/Tp, 
when  IITq  i  i  1/r^ 
when  f  >  l/f^. 


(4) 


We  refer  to  the  exponent  E  as  the  power  spectral  density  exponent.  Further,  JCF  established 
that  the  power  sp>ectral  density  exponent  is  determined  by  its  counterpart  in  the  distribution  of 
lifetimes,  viz.. 


1 


{0,  when  a  <  -3, 

-(3+ a),  when  -3  s  o  <  -1,  (5) 

-2,  when  a  ^  -1. 

Note  the  appearance  of  critical  frequencies  l/Tg  and  1/to  at  which  the  frequency  dependence  of 
the  PSD  changes  form.  Equations  (4)  and  (5)  constitute  the  JCF  connection  between  the 
distribution  of  lifetimes  and  the  PSD. 

The  JCF  connection  pertains  to  a  wider  range  of  systems  than  those  resulting  from  SOC 
processes.  The  results  are  consequences  of  the  absence  of  characteristic  length  and  time  scales 
and  therefore  apply  to  systems  that  exhibit  fractal  scaling,  etc.,  independent  of  an  organizing 
principle. 

Furthermore,  the  JCF  connection  does  not  depend  on  the  specific  form  of  the  cutoff 
power-law  distribution  of  lifetimes.  This  is  important  because  the  data  are  not  necessarily  well- 
described  by  exponentially  cutoff  power  laws.  For  example,  the  behavior  of  real  sandpiles  (ref  3) 
and  the  Barkhausen  effect  (ref  4)  indicate  that  the  distributions  are  cut  off,  but  the  Jest  form  is 
not  obvious.  The  distribution  of  lifetimes  in  the  Barkhausen  effect  in  three  ferromagnetic  metals 
(Metglas  2605S  (Fe7gB,jSi„),  poiyciystalline  iron  thermocouple  wire,  and  Alumel  (N^jAjMnj)) 
was  determined  to  be  better  represented  as  sharp  cutoff  than  exponentially  cutoff  weighted 
distributions  of  lifetimes  in  Reference  4. 

In  this  report  we  demonstrate  that  a  sharp  cutoff  weighted  distribution  of  lifetimes  gives 
rise  to  the  JCF  parameter  connections  with  the  cutoff  lifetime  (largest  lifetime  in  a  “size  effect" 
limited  distribution)  playing  the  role  of  Tg  in  the  exponentially  cutoff  distribution,  as  was  claimed 
in  Reference  4. 


THE  SHARP  CUTOFF  FORM 


Assume  that  the  weighted  distribution  of  lifetimes  C(T)  varies  (approximately)  as  a  sharp 
cutoff  power  law 


GO) 


0,  when  T  < 

'  r*,  when  Tq  i  T  i 
0,  when  T  > 


so  that  from  Eq.  (2),  the  PSD  is  given  by 


(«) 


2 


To 

S(f)  =vf  dr 
<0 


,|sm(it/r)J 


(7a) 


n/To 

=  v(Tty)-<“*5>  f  dx  x*(sm2(x))  (7b) 

=  v(ic/)-<-»>fK(a,ii>So,it/ro)  (7c) 

where  Eq.  (7c)  serves  to  define  the  function  W 

b 

m^Ab)  =  f  dx  x'sin^Cjc)  (7d) 


One  sees  that  S(f)  is  proportional  to  f^‘*^W(a,‘jtfio,TrfTo). 


ANALYTIC  RESULTS 


General  Solution  of  the  Sharp  Cutoff  JCF  Probiqn 


Mathematica  is  well-suited  to  the  evaluation  of  the  sharp  cutoff  expressions  for  the  PSD. 
For  general  values  of  a,  W{aMjb)  may  be  expressed  in  terms  of  generalized  hypergeometric 
functions  ^,(z).  After  some  editing,  the  Mathematica  input 


Integrate^  alpha5in{jc]  2, {x;a,bl]  yields; 


WidMfb) 


2(l+o)  {2  21  i2'2  2}  2(l+a)  I2  2}  \2*2  3/ 


,i*« 


a-fi 


It  is  clear  that  one  can  express  the  general  a  form  of  W{a4iJ>)  in  terms  of  functions  dependent 
on  the  limits  separately,  i.e., 

*  H»(«^)-w(a/i)  (*■) 


H<a,r)  = 


(») 


3 


Mathematica  2.1  evaluates  the  generalized  hypergeometric  function 

jc  <  20 

but  fails  for  larger  values  of  x.  (Mathematica  2.1  is  employed  to  obtain  numerical  values  of 
W{a^J))  for  arbitrary  a,  a,  and  b  via  NIntegrate{;e"alpha  S*n(»]*2,tefl,f»}].) 

Special  values  of  a 

When  a  is  integer  or  half-integer  valued,  W  may  be  expressed  in  terms  of  Fresnel 
integrals  and  cosine  and  sine  integrals,  a  better  known  set  of  hinctions,  which  were  recently  the 
subject  of  a  "Numerical  Recipes"  column  of  Press  and  Teukolsky  (ref  5). 

Employing  the  Mathematica  code: 

For[i  =  -8,  i  <  3, »  +  +, 

w[i/2,t_]:  =Evaluate[ 

MapfSimplify  [  Apart[#]]&,Integrate[t  *  (i/2)Sin[t]  *  2,t]]]]; 

one  obtains  the  following  expressions  w[a4\  for  integer  and  half-integer  a  values  in  the  range  of 
interest 


.r  ..  *1  _  -I  (I -2f^)co8(2f)  Stii<2r)  2SmihtegralCtt) 

i 


(9a) 


/  _\ 


wi-iM 


^  ^  ^  (3-16f^)co«(2f) 


32^Fresneld 


Ut 


j  _  4«ni(2/) 

15  1 

15f* 


(9b) 


w[-3^]  =  — ^  +  CosIfUegralOt)  - 


(9c) 


^-5M  = 


.  «>t(2r) 

1  3 

3f  ^  3f  ^ 


S^Fresne/cj 


yn 


4«ii(2f) 

3/^ 


(9d) 
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+  SinlmegraKZt) 


tK-2^1  =  *  ^5^ 

2t  2t 


yv[-3M  =  --1  +  +  2y/^  FresneJ^ 

/t  yft  \yfK, 


h{-1^]  =  -CosInUgralCh)  ^  log(r) 
2  2 


w[-l/2^1  = 


^  -  FresnelC 

( 

2it 

MOA 


=  i.  _  sm(2/) 


Mll/Vl  •  ^  t 
HU]  = 


i/n  Fresnels 

_ 


^sin(2r) 


8  4 

_  cos(2r)  rsin(2<) 

4  8  4 


(5e) 

(91) 

(9g) 

(9h) 


(9i) 


(id) 


(9k) 


With  these  analytic  expressions,  one  can  consider  various  ranges  of  /  and  demonstrate 
that  the  JCF  connections  obtain  for  the  special  a  values.  We  consider  a  few  cases  in  detail; 

A.  a  =  -2. 

1.  The  behavior  of  w[-24\  given  in  Eq.  (9e). 

a.  Small  t.  One  could  use  the  expansion  of  Eq.  (16)  of  Reference  5,  etc.,  to 
obtain  the  small  t  form  of  w[-24].  Here  we  employ  the  Series  command  from  Mathematica  to 
obtain 


h1-2j]  ^  — 

9  225 


2205 


+  (Xff 
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Large  t.  Using  the  large  r  limiting  form  (which  can  be  deduced  in 


b. 

Mathematica): 


Sinlnugral[t]  +  -y- 


-  it/2 


one  obtains: 


M.-2A  - 


It 

2 


2.  The  behavior  W[-2^J>]  and  the  PSD  exponent.  From  Eq.  (7c), 

a.  a  <<  b  <<  1  or  f  <<  I/Tq  <<  lAo- 

W[-2/iJ}]  «  w[-2,b]  -*  ic/Tq  and  SiJ)  -  const 

b.  a<<l<<bor  HT^  «f« 

W[-2/iM  «  H-2Jf]  -  nf2  and  S(f)  oc  /  > 


c.  a  <<  b  <<  1  or  l/Jo  <<  1/to  << /• 


The  7r/2  terms  in  H'[-2,a]  and  w[-2J>]  cancel  and 


1 

2Ttft^ 


Thus, 


s(/)  «/■*  X  nf  =  f  -^ 

Thus,  the  JCF  connections  are  established  for  the  case  a  =  -2. 

B.  a  =  -3/2. 

1.  The  behavior  of  w[-3/2,/]  given  in  Eq.  (9f). 

a.  Small  t.  One  could  use  the  expansion  of  Eq.  (10)  of  Reference  5,  etc.,  to 
obtain  the  small  t  form  of  w[-3/2,/].  Here  we  employ  the  Series  command  from  Mathematica  to 
obtain: 


6 


obtains: 


3 

2#  ^ 

Mi-mjt]  -  ^ 


1  11 
2t^  +  ^ 

21  495 


+  0(0" 


b.  Large  /.  Using  the  large  t  limiting  form  of  Eq.  (15)  of  Reference  5,  one 
M.-3M - -  *  yfH  ^  yfH 

i* 

2.  The  behavior  W[-3H^b]  and  the  PSD  exponent.  From  Eq.  (7c) 

a.  fl<<h<<lor/<<  l/To  <  <  l//o- 

W[-3I1^M  “  hI-3/2^]  -  |(it/ro)^  and  S{f)  -  const 

b.  a<<l<<i>or  VT^  «f«  l/r^. 

W[-3pi/iM  «  w(-3/2^1  -  yfn  and  S(f} 

c.  a<<h<<lor  I/Tq  <  <  1/to  <  <  /• 

The  \/ir  terms  in  w[-3/2^]  and  ^[-3/2^]  cancel  and 

WI-3PA*]  -  -L  -  ±  .=  1/y^ 
ya  ^ 

Thus, 

5(0  “  /  X 

Thus,  the  JCF  connections  are  established  for  the  case  a  =  -3/2. 
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C.  a  =  1. 

1.  The  behavior  of  ^[1^]  given  in  Eq.  (9k). 

a.  Small  t.  Employ  the  Series  command  from  Mathematica  to  obtain 

1 

MU]  -  0(0* 

8  4  18 

b.  Large  t. 

MUl  -  - 

4 

2.  The  behavior  W[\^Jb\  and  the  PSD  exponent.  From  Eq.  (7c), 

S(fi  «/-<••»»'(«, 

a.  a<<fc<<lor/<<  l/To  <<  lAo.  The  1/8  terms  in  w[\^]  and  w[\Jb] 

cancel  and 

W[\aM  »  (it/r//4  and  5(0  « /"*  X  /  -  const 

b.  1  <<  or  l/To  <</ (b>)th  subranges  included). 

W\\^M  «  Ml^l  -  and  S(f)  «  f-*  X  f  - 

Thus,  the  JCF  connections  are  established  for  the  case  a  =  1. 

Similar  analysis  can  be  applied  for  all  the  forms  in  Eq.  (9)  and  for  integer  and  naif- 
integer  a  generally. 
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NUMERICAL  RESULTS 


The  analytic  results  of  Eqs.  (7)  through  (9)  can  be  used  to  compute  S(‘nfTo)  for  specific 
values  of  a  and  Tg/to.  Figure  1  presents  typical  PSD  versus  frequency  results,  which  were 
obtained  for 


Tq  =  X  fp 

and  a  G  {-7/2,  -5/2,  -2,  -3/2, 1}.  The  value  of  a  increases  from  -7/2  for  the  top  curve  to  -t-1  for 
the  bottom  curve.  The  vertical  line  at  In('n^o)  =  12  corresponds  to 

In(ir^o)  =  -  12  -  12  -  12  =  0 

Thus  the  breaks  in  slope  occur  for 

ln(it^Q)  =  0  and  =  0 

as  advertised.  All  curves  become  /-independent  for  ln(ii/T o)  <  0  and  exhibit  inverse  square  PSD 
for  ln(T3;^o)  >  0-  The  top  curve  is  typical  of  a  <  -3  cases.  The  three  middle  curves  represent  the 
-3  <  a  <  -1  range  for  which  E  =  -(3-t-a);  the  variations  in  slope  are  apparent.  The  lowest  curve 
is  typical  of  the  a  >  -1  range.  Although  it  is  unlikely  that  one  could  observe  such  effects  in 
actual  PSD  curves,  the  "bumpiness"  to  the  right  of  the  high  frequency  transitions  is  real  (i.e.,  not 
numerical). 


ln(S(f)) 


Figure  1.  Power  spectral  density  versus  frequency  for  To  =  tg. 
Curves  are  shown  for  lifetime  distribution  exponent 
a  e  {-7/2,  -5/2,  -2,  -3/2, 1}.  a  increases  from  -7/2  for 
the  top  curve  to  +1  for  the  bottom  curve.  The  vertical 
line  at  ln(TifTo)  =  12  corresponds  to  ln('n/ro)  =  0. 
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Results  obtained  via  Nlntegrate  at  arbitraiy  a  tall  neatly  between  the  curves  shown  ir 
Figure  1  and  those  for  half-integer  and  integer  a  are  indistinguishable  from  curves  obtained  from 
the  analytic  expressions. 


CONCLUSIONS 

The  connection  between  the  distribution  of  lifetimes  and  the  PSD  established  by  Jensen, 
Christensen,  and  Fogedby  (ref  2)  for  the  case  of  exponentially  cutoff  distributions  of  lifetimes  has 
been  shown  to  apply,  with  natural  parameter  correspondences,  to  sharp  cutoff  distributions  of 
lifetimes.  Since  the  range  of  cutoff  forms  between  exponential  and  sharp  is  broad,  the  present 
results  suggest  that  the  JCF  cormections  will  obtain  to  a  very  wide  range  of  size-effect  modified, 
self-organized  critical  systems. 

The  power  of  symbolic  computational  ^sterns,  such  as  Mathematica,  is  nicely  illustrated 
by  the  present  analysis.  As  presently  described,  all  calculations  can  be  achieved  in  Mathematica. 
One  could  also  employ  the  analytic  results  as  starting  points  for  other  (e.g,,  Fortran  or  C) 
computer  programs. 
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